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Summary 

The  previous  work  of  the  author  in  which  the  suhcritical  flow  past  ellipsoids  was  calculated  is  extended  to 
supercritical  flows.  The  same  ellipsoidal  coordinate  system  is  used,  the  body-surface  boundary-condition  is 
then  applied  exactly  (in  the  numerical  sense).  By  means  of  a transformation  of  one  of  the  coordinates,  the 
infinite  flow  field  is  brought  into  a finite  space  for  the  calculation. 

The  complete  continuity  equation  is  approximated  by  the  usual  central  differencing  in  the  elliptic 
(suhcritical)  regions,  whilst  in  the  hyperbolic  (supercritical)  regions,  the  combination  of  central  and 
non-central  differencing  as  suggested  by  Albone  and  Jameson  is  used,  in  order  to  model  the  absence  of 
upstream  propagation  of  disturbances. 

Although  shock  waves  appear  in  the  calculations,  their  shape  and  position  is  only  approximately  determined 
las  e.g.  in  transonic  small  perturbation  theory)  since  the  difference  scheme  only  ensures  continuity  of  the 
potential  across  the  shock,  the  Rankine-Hugoniot  equations  not  being  satisfied.  A number  of  results  is 
presented,  for  flows  aligned  along  the  major  and  second  axes,  and  also  yawed  relative  to  these  two  axes. 

The  work  reported  here  was  done  in  connection  with  the  link  between  the  R.A.E.  and  the  University  of 
Southampton. 


* Now  at  the  Department  of  Mathematics,  Imperial  College, 
t Replaces  R.A.E.  Technical  Report  750X7— A R C.  36  304. 
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1.  Introduction 


The  development  of  powerful  electronic  computers  in  the  lust  few  years  has  enabled  numerical  techniques 
to  he  developed  to  deal  with  the  problem  of  transonic  flow  past  two-dimensional  aerofoil  sections,  and  more 
recently  past  three-dimensional  bodies. 

Magnus  and  Vosh'hara 1 overcame  the  difficulty  of  the  presence  of  both  suhcritical  and  supercritical  regions 
within  the  flow,  by  considering  the  steady-state  solution  to  be  the  limit  of  a non-steady  solution,  and  thereby 
ensured  a fully  hyperbolic  system,  at  the  expense  of  introducing  into  the  problem  a further  dimension,  namely 
time. 

Nlurman  and  Cole~  introduced  the  mixed  finite  difference  scheme.  This  approximates  derivatives  in  elliptic 
regions  in  the  usual  way,  by  central  differences,  whilst  in  hyperbolic  regions,  derivatives  in  one  of  the 
coordinate  directions  (the  one  most  closely  aligned  with  the  flow)  are  approximated  using  backward 
differencing  (with  respect  to  the  flow  direction).  Other  derivatives  are  centrally  differenced.  Physically  this 
scheme  is  an  attempt  to  model  the  absence  of  upstream  propagation  of  the  flow  information  in  the  supersonic 
zones,  disturbances  at  a point  being  transmitted  solely  within  the  downstream  Mach  cone  w hich  has  its  apex  at 
that  point. 

This  method  was  also  applied  successfully  by  Garabedian  and  Korn1  for  two-dimensional  aerofoils. 
However  if  the  local  velocity  in  the  supersonic  region  departs  appreciably  from  one  of  the  coordinate 
directions,  numerical  instability  can  occur.  To  overcome  this  problem,  Albone4  and  Jameson'  have  suggested 
transforming  the  principal  part  of  the  potential  equation,  in  hyperbolic  regions,  into  a new  coordinate  system, 
with  one  axis  aligned  along  the  local  flow  direction.  Second  derivatives  in  the  direction  of  the  flow  are  then 
differenced  backwardly,  whilst  other  derivatives  are  centrally  differenced.  Elliptic  regions  are  treated  in  the 
usual  way. 

This  report  is  an  extension  of  the  work  on  suhcritical  flow  past  ellipsoids  of  Duck'’.  The  flow  around 
ellipsoids  without  circulation  is  considered  (implying  no  trailing  vortex  sheets)  using  the  full  equations  of 
motion  and  the  ellipsoidal  coordinate  system,  which  facilitates  the  exact  implementation  of  the  body-surface 
boundary-condition.  Further,  by  means  of  a simple  transformation,  the  infinite  physical  flow  field  is  brought 
into  a finite  working  space.  In  the  supersonic  regions,  the  proportion  of  retarded  differencing  suggested  in 
Refs.  4 and  5 is  used. 

Attention  is  restricted  to  local  Mach  numbers  only  slightly  in  excess  of  unity  since  a velocity  potential  is  used, 
ignoring  any  changes  in  entropy,  and  any  rotation  in  the  flow. 

2.  Equations  of  Motion 

Throughout  we  shall  use  the  ellipsoidal  coordinate  system  (f ',  £3)  where 

jr  =[H-(^1)2]4  dn  f 2 sn  g\ 


y = [cr  + (f 1 )2]J  cn  ^ cn  f ' 

and 

2 = f 1 sn  dn 

where,  for  conciseness,  we  write 

sn  =sn  s/a-) 


and 

sn  £'  = sn  (f  \ V 1 - a) 

(and  similarly  for  the  other  Jacobian  elliptic  functions), 

rr  is  a parameter  that  partly  determines  the  ratio  of  axes  of  the  ellipsoid,  such  that 


0 < a < 1 . 


The  continuity  equation  in  the  transformed  ellipsoidal  system  is  now 

i r.  i.  i r.  <t>2  1. 
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Further,  in  order  to  determine  the  speed  of  sound  a,  we  use  the  Bernoulli  energy  equation 

a1  . 1 , 


...  a2«t> 

" a(f)2 


!'./  = 2,3. 


The  flow  speed  is  expressed  by 


- + -q~  = constant. 

y- 1 2 


2 4>i  ^ 4>;  , 4>5 

q=tf*A?A\ 

= u\  + u\  + mT 


To  complete  our  description  of  the  ellipsoid,  on  which  the  boundary  condition  is  satisfied  we  set 

£‘=£0  (a  constant) 

= Vcrtn  £0 


x~  i y z~ 

1 +(f,')2+rr  + (^)5  + (iI) 


On  this  surface  we  apply  the  body-surface  boundary-condition 

4>,  =0. 

Since  we  are  concerned,  throughout,  with  bodies  without  circulation,  the  boundary  condition  at  infinity  is 

4>  = 4>r  at  £' = K\i' =cc), 
here  4>T  is  the  undisturbed  (freestream)  velocity  potential. 
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In  this  coordinate  system,  t.long  the  hyperbolae  £~  = ±K,  £'  = ±K‘,  we  have  B\  — 0,  implying  a vanishing 
Jacobian  of  the  transformation.  It  was  shown  by  Duck*’  that  along  these  lines  the  potential  equation  reduces 
simply  to 

+<!>,,  = 0,  (2) 


and  the  flow  velocity  is  given  by 


<1>I 

021  ^ 3 

_l_  *• 

[or  + (£')2; 

]’  rr(  1 -tr)[tr  + (£')'] 

4>f 

tfi2.,  + <l>2, 

[<7  + (£')2'. 

]’  (T(l  -cr)[(r  + (£'):) 

A fuller  description  of  the  ellipsoidal  coordinate  system  is  given  in  Ref.  6. 


3.  Numerical  Techniques 


3.1.  Finite  Differences  Equation 

We  follow  the  example  of  Albone4  and  Jameson"  by  rearranging  the  principal  part  (that  is  all  terms  involving 
second  derivatives)  of  the  potential  equation  ( 1 ) into  canonical  form,  in  the  following  way: 


( 1 -M2)<l>„+V2<l>-<t>„ 


where 


<t>„  = — 
<T 


ul  . ul  ^ 2u \U z . , 2 u2u, 

A I A%  Ay  A,A:  AtAy  A2Ay 


V <!>  = i 


A; 


A i 


and  M = q/a  is  the  local  Mach  number. 

The  first  order  derivatives  of  1 1 ) remain  unchanged. 

For  the  purposes  of  computation,  the  total  potential  <t>  is  now  split  up  in  the  form 


(J>  = <J)X  + <t>, 


where  <t>,  denotes  the  freestream  (undisturbed)  potential,  and  is  a known  quantity,  with  a singularity  at 
infinity — we  solve  only  for  <t>.  the  perturbation  (unknown)  potential. 

The  partial  differential  equation  is  now  ready  to  be  differenced. 

Assuming  lm  + 1,  m„  + 1.  n„  + 1 points  in  the  £\  £:  and  £ 1 directions  respectively,  then  the  cell  dimensions 
in  the  differenced  system  (for  the  general  case),  are  given  by 

\£'=(K'-£?,) 

A£:  = 4 K/m„ 

and 

\£x  = 2K'/nm. 

The  coefficients  of  K and  A,"  in  the  above  expressions  for  \£  and  ±£ ’ may  be  reduced  if  any  symmetry  is 
present  in  the  problem. 
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All  first  derivatives  of  <t>  are  differenced  centrally  throughout  the  entire  flow  field,  for  example 


n 


* t (fit.rn.n  1 

2Af'~ 


+ 0((A£')2) 


(where  /,  m,  n identify  the  point  in  the  field). 
Second  derivatives  of  </>  contained  within  V<t>- 


ct>„  are  also  centrally  differenced,  for  example 


<£/. 


'~(j> 

afr 


+ 0((MY) 


and 


(•^nj/.m.n 


<t>l.m  + l.n*l  ~ <£/.m  -l.n  + 1 ~ <t> 


l,m  + 1 ,n  1 + <£  l.m-l.n  1 


4Af2Af2 


+ 0((A^:);(A^')2). 


For  elliptic  (subcritical)  points  the  second  derivatives  of  </>  contained  within  (1  are  differenced 

centrally,  as  above.  However  for  hyperbolic  (supercritical)  points,  the  differencing  is  non-central,  for  example 


tri.n  ± 1 + </>, 

,m,nt  2 , 

— . . *1.2  •"  C/(.l£  ), 

(Af  ) 


where  the  positive  (negative)  sign  is  taken  for  <t>,  negative  (positive). 
Similarly 


(4>2i)l.m.n  =sign  (<t>2*t>v) 


l,m,n  t.  1 


+ 4>i. 


if  if 


+ CHAf  + if) 


where  the  signs  are  chosen  depending  on  the  signs  of  <t>:  and  <t>i. 

Along  the  singular  lines,  we  use  equation  (2),  which  is  elliptic,  for  both  subcritical  and  supercritical  points, 
central  differencing  always  being  applied  to  this  equation.  Points  on  either  side  of  these  singular  points,  in  the 
f and  f directions  are  in  fact  equivalent,  due  to  the  periodicity  of  the  working  space  (see  Fig.  1). 

Symmetry  arguments  may  be  invoked  (if  present)  in  order  to  reduce  both  computer  time  and  store 
requirements.  For  flows  with  freestream  velocities  aligned  along  one  of  the  axes  of  the  ellipsoid,  two  quadrants 
ire  required  for  supercritical  calculations  (the  fore  and  aft  anti-symmetry  of  subcritical  cases  is  no  longer 
present). 

Referring  to  Fig.  1,  if  <t>x  =x,  any  adjacent  combination  of  quadrants  A and  B (or  C and  D)  may  be  used, 
symmetry  may  then  be  invoked  along  the  perimeter  of  each  plane  f = constant,  because  of  the  up  and  down, 
ind  side  to  side  symmetry  present.  If  <t>X)  = y,  any  adjacent  combination  of  quadrants  A and  C (or  B and  D)  may 
ised.  In  addition  to  symmetry  considerations,  the  periodicity  of  the  working  space  must  be  taken  into  account. 
As  an  example,  consider  the  plane  = constant  over  the  region 

O^f^K 


0«f’ SS2AT' 

is  shown  hatched  in  Fig.  1 . Symmetry  may  be  invoked  along  OQRSO,  but  not  along  OTO.  But  because  of 
seriodicity,  a point  such  as  V , described  in  ellipsoidal  coordinates  by  if , K + (')  maps  to  the  same  point 

n cartesian  (physical)  space  as  the  point  (f K-if,  2K'-f)  in  the  working  space.  Similar  arguments  apply 
o other  combinations  of  B and  D (or  A and  C). 

For  yawed  cases,  more  quadrants  must  be  brought  into  the  calculations,  with  arguments  similar  to  those 
;iven  above  invoked  in  order  to  ‘close’  the  working  space. 

Finally  we  apply  the  body  surface  boundary  condition  by  means  of  reflection. 

If  H,  1 « K'  as  1 /«  lm  + 1 then 

m.n  ^ 2 Af  (4^x  1 Ifsodv  surface. 
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This  then  gives  values  of  </>  for  the  set  of  points  situated  just  inside  the  body  that  are  required  in  the  difference 
equation. 

3.2.  Solution  of  Difference  Kquation 

The  difference  equation  may  now  be  expressed  in  a form  suitable  for  line  relaxation.  For  solution  along  lines 
£',£'  = constant  (spokes),  the  difference  equation  may  be  expressed  in  standard  tridiagonal  form 

1 .m.n  = «/  for  1 =£ /«/„,(</>,„,»  j 0) 


where  b,nt  = 0 and  </,  = 0. 

As  many  terms  as  possible  are  placed  on  tjie  left-hand  side  of  the  above  equation,  the  remainder,  including 
all  the  second  order  derivatives  of  the  freestream  contribution  to  the  potential  and  also  most  of  the 
contributions  in  the  difference  equation  from  other  mesh  points,  being  placed  on  the  right-hand  side. 

The  resulting  matrix  equation  for  <t>  may  be  solved  using  the  well  known  Choleski  method  (see  for  example 
Varga  ). 

Having  solved  for  <b  in  this  way,  we  apply  relaxation  to  update  the  old  value  of  d> 

(b'in, + ( 1 ~ 

being  obtained  from  the  block  inversion  described  above,  and  superscript  (.V)  denotes  the  value 
obtained  (after  relaxation)  from  the  Nth  iteration. 

The  relaxation  parameter,  w,  was  usually  set  at  1 -8  in  elliptic  regions,  and  0-96  in  hyperbolic  regions  at  the 
start  of  the  computation,  and  was  reduced  if  instability  occurred.  Frequently  the  computation  would  undergo  a 
period  of  instability,  during  which  the  relaxation  parameters  were  reduced.  The  numerical  process  would  then 
become  stable  again,  and  the  relaxation  parameters  could  then  be  increased  without  further  instability . 

The  direction  of  ‘sweeping-  of  the  lines  along  which  the  solution  was  obtained  was  varied,  including 
following  the  flow  and  directly  opposed  to  the  flow,  with  identical  results. 

The  region  close  to  the  shock  was  usually  very  slow  to  convergence  during  the  iteration  process.  Often  the 
solution  in  this  region  would  still  be  changing  markedly,  whilst  in  the  rest  of  the  flow  field,  the  solution  had 
settleu  down.  For  this  reason  convergence  was  usually  assumed  when  the  maximum  change  in  potential, 
anywhere  within  the  field  was  less  than  1x10  h per  iteration.  This  implies  convergence  of  Mach  number  of 
approximately  5x10  * per  iteration. 

The  shock  region  also  appears  the  most  sensitive  to  changes  in  mesh  size  as  shown  in  the  example  in  Table  1 
which  refers  to  the  case  illustrated  in  Figs.  2 to  6 — presumably  this  is  a result  of  the  large  gradients  of  the 
potential  around  the  shock  (situated  near  y =0-45).  Table  1 indicates  that  accuracy  requires  a particularly 
large  number  of  points  in  the  £'  direction  which  approximately  corresponds  to  the  direction  normal  to  the 
shock.  Most  calculations  for  the  unyawed  cases  were  performed  using  a 2 1 x 2 1 x 41  mesh  in  the  £ ',  £:  and  £ 
directions  respectively,  whilst  those  for  yawed  cases  were  made  on  a 2 1 x 4 1 x 4 1 mesh.  With  these  two  grids, 
the  Mach  numbers  appear  in  general  accurate  to  within  about  l per  cent. 

To  reduce  computation  time,  the  calculation  was  started  on  a coarse  grid  and  continued  on  a medium  grid  (of 
half  the  original  cell  dimensions),  before  being  transferred  to  a fine  grid  (of  a quarter  the  original  cell 
dimensions).  Supercritical  examples  take  about  10  times  the  number  of  iterations,  and  about  20  times  the 
computer  time  of  comparable  subcritical  examples.  A typical  unyawed,  supercritical  case  takes  about  19 
minutes  on  a CDC  7600  computer,  with  100  iterations  on  the  two  coarse  grids,  and  600  on  the  fine  grid.  Yawed 
supercritical  cases  take  about  double  the  number  of  iterations  and  four  times  the  computer  time  of  comparable 
unyawed  examples. 


4.  Results 

By  means  of  a number  of  linear  interpolation  routines,  the  converged  solutions  were  used  to  produce 
various  Mach  number  distributions  in  cartesian  (jc,  y,  z)  space,  both  on  the  body,  and  also  out  into  the  flow 
field.  Several  of  these  distributions  are  presented  in  Figs.  2 to  33  (as  shown  listed  in  Table  2).  It  should  be 
realised  that,  because  of  the  nature  of  the  interpolation,  the  distributions  away  from  the  body,  in  particular, 
should  be  treated  with  a certain  amount  of  caution,  because  of  the  lower  density  of  points  as  infinity  is 
approached  in  the  physical  field. 

Figs.  2 to  1 1 show  distributions  of  Mach  number  around  two  different  bodies,  with  flows  aligned  along  the 
second  major  axis  (<t>,  =y).  The  first  body  (hereafter  referred  to  as  body  I)  has  a ratio  of  axes 


1 : 0-7089: 0-0709  in  the  a : v:z  directions  respectively,  whilst  the  second  body  (referred  to  as  body  II)  has  a 
corresponding  ratio  of  axes  of  1 : 0-2010  : 0-0200.  Figs.  12  to  16  show  the  flow  around  body  I with  a free- 
stream  flow  aligned  along  the  major  axis  (<t> *,  = a).  For  all  these  examples,  it  appears  that  the  Mach  number 
across  the  maximum  'span'  of  the  body  is  approximately  constant,  although  the  examples  of  higher  aspect  ratio 
exhibit  a slight  increase  in  Mach  number  in  the  region  of  the  tips  ( see  Figs.  4.  9,  14).  The  exact  incompressible 
solution  predicts  constant  flow  speed  across  the  maximum  'span'  of  the  ellipsoid.  For  the  high-aspect-ratio 
body  (Fig.  8)  there  appears  to  be  a tendency  for  the  position  of  the  maximum  Mach  number  to  move  away  from 
the  plane  of  symmetry  and  consequently  there  may  be  an  increase  in  shock  strength  as  the  'tips'  are 
approached.  For  the  low-aspect-ratio  body  (Fig.  1 3)  the  opposite  occurs  and  away  from  the  plane  of  symmetry 
the  shock  wave  appears  to  disperse  progressively  into  a distributed  compression.  Figs.  6,  11,  16  show  the  flow- 
field  in  depth  in  the  plane  of  symmetry  and  emphasize  the  increasing  departure  of  the  flow  from  fore-and-aft 
symmetry  as  the  freestream  Mach  number  increases.  Further  the  shocks  appear  to  meet  the  body  surface 
approximately  perpendicularly.  Figs.  3,  13  indicate,  in  the  strongest  regions  of  the  shock,  after  an  initial  rapid 
compression,  a very  weak  expansion,  followed  by  a second  compression.  (This  effect  is  most  marked  in  the 
example  of  Fig.  13.)  For  reference,  the  positions  of  the  grid  points  along  the  centreline  are  shown  in  these 
figures.  All  other  results  were  obtained  by  interpolation.  Along  the  centre  line  the  shock  is  smeared  over  two 
or  three  mesh  points.  It  is  expected  that,  because  of  the  more  favourable  distribution  of  'calculating'  points  for 
the  <1*1  = a example  (the  shock  occurs  in  the  region  of  one  of  the  singular  points  on  the  body,  a region  in  which 
there  is  a natural  'bunching'  of  mesh  points  in  the  a -direction)  the  shock  definition  is  rather  better  than  in  the 
other  examples. 

Downstream  of  the  shock,  the  Mach  number  distributions  appear  to  almost  regain'  the  symmetry  of 
subcritica!  flows,  the  Mach  number  at  a position  behind  and  away  from  the  shock  matching  closely  with  that  at 
the  corresponding  position  upstream  of  the  shock.  Details  of  subcritical  flows  past  unyawed  ellipsoids  are 
given  in  Ref.  6. 

Distribution  of  Mach  numbers  for  flows  around  body  f , with  the  freestream  flow  inclined  at  45  degrees  to  the 
a and  y axes  (<t>x  = (x  + y)/s/2)  are  shown  in  Figs.  17  to  26.  Figs.  17  to  19  are  for  a subcritical  example 
(Mx  =0-8),  for  comparison  with  a supercritical  case  (M*  = 0-945),  details  of  which  are  giv  en  in  Figs.  20  to  26. 

The  structure  of  the  shock  for  this  supercritical  example  is  similar  to  that  of  the  unyawed  examples  discussed 
previously. 

The  acceleration  following  the  rapid  compression  this  time  results  in  a second  (just)  supersonic  zone,  before 
the  second  compression  (see  Fig.  21). 

The  position  of  the  stagnation  points,  in  comparison  with  the  subcritical  case  appear  little  changed,  in  spite 
of  the  appearance  of  a shock.  Again,  well  away  from  the  shock,  the  flow  appears  to  regain  the  symmetry  of 
Mach  number  distribution  present  in  the  subcritical  example  (Figs.  17  and  20). 

Figs.  27  to  33  refer  to  the  flow  around  body  II  with  freestream  flow  inclined  at  45  degrees  to  the  x and  v axes, 
with  Mt  = 0-97. 

The  main  characteristic  of  these  figures  is  the  absence  of  any  noticeable  shock,  the  flow  appearing  still  to 
have  the  symmetry  associated  with  totally  subcritical  flows. 

Finally,  it  should  be  stressed  that,  since  the  Rankine-Hugoniot  relations  have  not  been  satisfied  across  any 
shock  (only  continuity  of  potential  is  satisfied)  then  the  actual  structure  and  position  of  any  shocks  arising  from 
these  calculations  should  be  treated  with  caution. 


5.  Conclusions 

It  has  been  shown  that  numerical  solutions  may  be  obtained  for  the  supercritical  flow  around  ellipsoids 
without  circulation  and  under  steady,  inviscid  conditions. 

Although  the  Rankine-Hugoniot  equations  have  not  been  applied,  the  full  equations  of  motion  (in 
ellipsoidal  coordinate  form)  have  been  used  and  the  boundary  condition  on  the  surface  of  the  ellipsoid  has 
been  applied  as  exactly  as  the  differencing  schemes  allow.  Using  a transformation  of  one  of  the  coordinates, 
the  entire  flow  field  has  been  brought  into  a finite  computation  space.  Shocks  appear  in  the  converged 
solutions,  generally  smeared  over  two  or  three  grid  points. 

The  numerical  scheme,  using  the  proportion  of  central  and  non-central  differencing  as  suggested  by 
Albone4  and  Jameson'  has  generally  been  reliable  with  only  intermittent  outbreaks  of  instability.  Computa- 
tion times,  however,  are  still  appreciably  longer  than  their  subcritical  counterparts,  primarily  because  of  the 
slow  convergence  at  points  in  the  neighbourhood  of  the  shock. 
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1 1ST  OF  SYMBOLS 


a 

iii,  bi,  Ci,  di 

A, 

B, 

K(K') 

/pit  ' 1 ’ 

M 

<■ I 

U, 

x,  y, z 
y 

lUTf1 

f 

(T 

<i> 

«b 

<t>x 

<t>, 

*« 

4>„ 

Subscripts 

i,  /,  * 
/,  w,  n 
0 

oo 

Superscripts 

i,i,k 

N 


Speed  of  sound 

Coefficients  in  matrix  equation 
Curvilinear  metrics 
Magnitudes  relating  to  the  A, 

(Complementary)  complete  elliptic  integrals  of  the  first  kind 
Number  of  points  in  £ directions 

Mach  number 
Flow  speed 

velocity  component  in  g direction 
Cartesian  coordinates 
Ratio  of  specific  heats 
Ellipsoidal  coordinates 
Transformed  £ 1 coordinate 
Body  description  parameter 
Perturbation  potential 
Total  potential  +<f> 

Freestream  (undisturbed)  potential 

a<t>/af  (/  = 2. 3);  ao/a#1  (<  = l) 
a>/a^  af*  a = 2, 3);  a-<j »/af  af'  (<  = l) 

Second  derivative  in  streamwise  direction 

Indices  running  from  1 to  3 
Indices  identifying  points  in  field 
Body  surface  conditions 
Freestream  conditions 

Indices  running  from  1 to  3 
Value  obtained  from  Mh  iteration 
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TABLE;  1 


Variation  of  Mach  Number  on  Body  Surface  with  Mesh  Size  Along 
x = o,  «J>r  = y,  Mx  = 0-92,  tr  = 0-5 
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15x21  x41 


21  x 15x41 
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1-2014 
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TABLE  2 

Details  of  Calculations  Presented 


Body  I:  o-  = 0-5  £,’=0-1 

Ratio  of  axes:  1:0-7089:0-0707 

Body  II:  o-  = 0-04  £,  = 0-1 

Ratio  of  axes:  1:0-2010:0-0200 
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Figures 
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Freestream  conditions 
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Mx  = 0-92 
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4>x  = y 

Mr  = 0-88 

12-16 
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= X 

Mx  = 0-97 

17-19 
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<t>x=(x+y)/J2 

Mr  = 0-80 

20-26 

I 

4>r  = (x+y)/s/2 

Mr  = 0-945 

27-33 

II 

4>x  = (x+y)/v/2 

Mr  = 0-97 
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Fk..  2.  I ines  of  constant  Mach  number  on  body  surface  A/,  =0-924>,  =y  Bodyl. 


Mach  number  distributions  along  lines  y = constant  on  body  surface  (y  =£  0)  Mx  — 0-92  <J>X  = y Body  I. 


Fig.  5.  Mach  number  distributions  along  lines  y = constant  on  body  surface  (y  s* 0 ) Mx  = 0-92  «t>r  = y Body  I 


Fig.  6.  Lines  of  constant  Mach  number  in  the  plane  x = 0 Mx  = 0-92  = y Body  I. 


Fig.  7.  Lines  of  constant  Mach  number  on  body  surface  ,\fx  = 0-88<l\  = y Body  II. 


Fig.  8.  Mach  number  distributions  along  lines  x = constant  on  body  surface  Mx  = 0-88  <J>,  = y Body  II 


Fig.  10.  Mach  number  distributions  along  lines  y = constant  on  body  surface  (v  ^0)MX  =0-88  <t>»  = y Body  II. 


Mesh  points 


Fig.  13.  Mach  number  distributions  along  lines  y = constant  on  body  surface  M 


80 


Lines  of  constant  Mach  number  on  body  surface  M , = 0-80  <b  v 


Fi<;.  18.  Mach  number  distributions  along  lines  y = constant  on  body  surface  Mx  = 0-80  d>,  = (1/v  2)(x  + y )Body  I 


Fie..  19.  Mach  number  distributions  along  lines  .r  = constant  on  body  surface  M,  = 0-80  <t>x  = (1  /v  2)(.v '+  y ) Body  I 


Fit;.  20.  Lines  of  constant  Mach  number  on  body  surface  = 0-945<t> 


Fig.  21.  Mach  number  distributions  along  lines  v = constant  on  body  surface  (v  « 0)  M x = 0-945  <J>,  = ( 1/s  2 )(a-  + y ) Body  I. 


Fig.  22.  Mach  number  distributions  along  lines  y = constant  on  body  surface  (>•  >0 )MX  = 0-945  <£v  = (l/\  2)(x  + v)Body  I 


Mesh  points 


Fig.  23.  Mach  number  distributions  along  lines  x = constant  on  body  surface  (jr  s 0)  Mx  — 0-945  — (1/s  2)(x  + y ) Body  I. 


Fig.  24.  Mach  number  distributions  along  lines  x = constant  on  body  surface  (x>0)Mx-  0-945  <t>.T  = (1/V2)(x  + y ) Body  I. 


) Body  I. 


. 28.  Mach  number  distributions  along  lines  y = constant  on  body  surface  (>'  «0)M,  = 0-97  <t>  - (l/v'2)(.v  + v)  Body  II. 


Fig.  29.  Mach  number  distributions  along  lines  y = constant  on  body  surface  (y  >0) M,  = 0-97  <t>,  = (!/'  2)(.v  +y)Body  II. 


Fio.  30.  Mach  number  distributions  along  lines  x = constant  on  body  surface  (r  =£  0)  Mx  = 0-974)*  = (1  /v  2)(.r  + v ) Body  II 


Fig.  31 . Mach  number  distributions  along  lines  * = constant  on  body  surface  (jc  >0)  Mx  = 0-97  <J>*  = (1/v  2)(x  + y )Body  II. 


Fig.  32.  Lines  of  constant  Mach  number  in  the  plane  y =0  Mx  = 0-97  <|>T  = (1/v  2)(x  + y (Body  II. 
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y 


Fic..  33.  Lines  of  constant  Mach  number  in  the  plane  x = 0 Mx  = 0-97  <!>,  = (l/v/2)(jr  + y ) Body  II. 
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hyperbolic  (supercritical)  regions,  the  combination  of  central  and  hyperbolic  (supercritical)  regions,  the  combination  of  central  and 

non-central  differencing  as  suggested  by  Albone  and  Jameson  is  used.  non-central  differencing  as  suggested  by  Albone  and  Jameson  is  used, 

in  order  to  model  the  absence  of  upstream  propagation  of  distur-  in  order  to  model  the  absence  of  upstream  propagation  of  distur- 


